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ABSTRACT 



We present Herschel-PACS observations of rest-frame mid-infrared and far-infrared spectral line emissions from two lensed, ultra- 
luminous infrared galaxies at high redshift: MIPS J142824.0+352619 (MIPS J1428), a starburst-dominated system at z = 1.3, and 
IRAS F10214+4724 (F10214), a source at z = 2.3 hosting both star-formation and a luminous AGN. We have detected [O I]63yum and 
[OIII]52yum in MIPS J1428, and tentatively [Om]52fim in F10214. Together with the recent ZEUS-CSO [CII]158yum detection in 
MIPS J1428 we can for the first time combine [01], [CII] and far-IR (FIR) continuum measurements for photo-dissociation (PDR) 
modeling of an ultra-luminous (L;^ >10'^Lq) star forming galaxy at the peak epoch of cosmic star formation. We find that MIPS 
J1428, contrary to average local ULIRGs, does not show a deficit in [O I] relative to FIR. The combination of far-UV flux Go and gas 
density n (derived from the PDR models), as well as the star formation efficiency (derived from CO and FIR) is similar to normal or 
starburst galaxies, despite the high infrared luminosity of this system. In contrast, F10214 has stringent upper limits on [OIV] and 
[S III], and an [O III]/FIR ratio at least an order of magnitude lower than local starbursts or AGN, similar to local ULIRGs. 
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ISO and Spitzer observations have shown that the fraction of 
cosmic star formation in infrared luminous and ultra-luminous 
galaxies ( LIRGs/ULIRGs: L,fl=10"Lo to a few times IO'^Lq) 
increases from a few percent at z~0 to more than 50% at z>l 
(Genzel & Cesarsky 2000, Soifer et al. 2008). Herschel offers the 
first opportunity to study the population of high-redshift infrared 
bright galaxies at wavelengths where they emit most strongly, 
both through photometry and high-resolution spectroscopy. In an 
attempt to probe for evolution of the ISM in the infrared galaxy 
population and to explore the limits of the PACS spectrometer 
we are investing several tens of hours for extremely deep high- 
resolution FIR spectroscopy of a (small) sample of galaxies with 
L/j? >1O'^L0 at redshifts > 1. These data will be complemented 
by and compared to the results from our ongoing Herschel spec- 
troscopy of nearby ULIRGs, metal-rich and metal-poor galaxies, 
and starburst and AGN template objects. We will thus explore 
a wide range of parameter space and redshift. In this letter we 
present the first results on two high redshift objects, demonstrat- 
ing the capabilities of the PACS spectrometer in long integra- 
tions and high-z spectroscopic observations. 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 

** Visiting scientist at Max-Planck-Institute for Extraterrestrial 
Physics (MPE), Garching, Germany 



2. Targets 

We present spectra of two high-redshift galaxies (MIPS 
J142824.0H-352619 and IRAS F10214-I-4724) that were ob- 
served in our initial Science Demonstration Phase observations 
from two PACS GT Key Projects: SHINING (PI E. Sturm) and 
the "Dusty Young Universe" (PI K. Meisenheimer). Both galax- 
ies are gravitationally lensed bringing them potentially within 
the sensitivity of the PACS spectrometer in long integrations. 

IRAS F10214+4724 (hereafter F10214) is a ULIRG at red- 
shift z=2.2855. It is known to harbor a coeval AGN, with 
a starbursting host galaxy. Differential magnification of these 
components (//(AGN)///(host)=3, Eisenhardt et al. 1996) poten- 
tially complicates the interpretation of any emission detected. 
Applying a bolometric magnification factor of 12 (Ao et al. 
2008) for the host galaxy yields an intrinsic luminosity of 
Lf/fi ~7 X 10'^ Lq. Teplitz et al. (2006) presented a mid-infrared 
ultra-deep spectrum with Spitzer IRS at low spectral resolution. 
Although classified as a Type II AGN, the mid-IR spectrum 
shows strong silicate emission around 10/vm similar to Type I 
QSOs, and very weak (if any) PAH dust features traditionally 
used as tracers for star formation. Nevertheless, taking differen- 
tial extinction into account, the AGN might contribute less than 
50% to the bolometric luminosity of this object. 

MIPS J142824.0+352619 (MIPS J1428 hereafter) was dis- 
covered in Spitzer-MIY'S images of the NOAO Deep Wide Field 
Survey Bootes field. Subsequent Spitzer-lRSi low resolution 
spectra (Desai et al. 2006), sub-mm (Borys et al. 2006, lono et 
al. 2006, Hailey-Dunsheath et al. 2010) and VLA observations 
(Higdon et al. 2005) revealed that it is an extremely luminous. 



E. Sturm et al.: Herschel-PACS spectroscopy of IR-bright galaxies at high redshift 



rest wovelenqth \u.m] 
63.0 63.2 63.4 




3.5X10-" 



E 3.0x10"" 

E ,7 

■^- 2.5x10" 



•S 2.0x10"" 

0) 
T3 

J 1.5x10"" 
l.OxlO"" 



146.0 146.5 147.0 147.5 148.0 

51.6 51.8 52.0 52.2 



MIPS J1428 |— 1 
[0 III] Sl.BMm 






(b) 



120.0 



120.5 



121.0 



121.5 



51.6 






51. 


7 


51 


.8 




51 


.9 


52.0 


5.2x10"" 


■ F10214 












(c) ' 


5.0x10"" 


. [° 


III 


51 


.8/im 




















1 


"-\ 










4.8x10"" 










/ 




\ 






. 




rL 






rUi 




k 








4.6x10"" 








_ _ 


1 / 






iJ 


-- 


n 


4.4x10"" 


L 


j_ 


r-' 


. 










Ml 


4.2x10"" 






Lr 












■ 


4.0x10"" 


, 




, 


, 


_ 






, 


, 



1 69.6 1 69.8 1 70.0 1 70.2 1 70.4 1 70.6 1 70.8 1 71 .0 
observed wovelength [/im] 

Fig.l. Observed PACS spectra. Top: [O I]63.2;um in MIPS 
J1428. Overplotted is a Gaussian profile on top of a linear con- 
tinuum fit. Also shown (upper dotted line) is a NH3 model (ar- 
bitrarily scaled), and the residual (lower dotted line) after sub- 
tracting the line+continuumH-NH3 fit from the data. Middle: [O 
III]51.8//m in MIPS J1428. Overplotted is a gauss profile on 
top of a linear continuum fit. Bottom: [O III]51.8//m in F10214. 
Overplotted is a gauss-i-linear continuum fit. 



starburst-dominated galaxy at z=1.325. This galaxy shows no 
signs of an unobscured or obscured AGN. It is magnified by a 
foreground galaxy, with a relatively small magnification factor 
of;U<8(Ionoetal. 2006). 



3. Observations and Data Reduction 

The observations were taken with the PACS spectrometer 
(Poglitsch et al. 2010) on board the Herschel Space Observatory 
(Pilbratt et al. 2010) in high resolution range spectroscopy mode. 
Most of the data reduction was done using the standard PACS 
reduction and calibration pipeline (ipipe) included in HIPE 2.0 



I340Q, together with some additional steps to correct for imper- 
fect offset compensation by the chop/nod observing technique, 
which can be important at these faint signal levels. The contin- 
uum in each of the 16 spectral pixels was scaled to the median 
value to correct for residual flat field effects. Finally the two nod 
positions were combined to completely remove the sky (tele- 
scope) background. Given that our targets are point sources to 
PACS we have measured line fluxes from the spectrum in the 
central spatial pixel of the 5x5 pixel FOV of the PACS spec- 
trometer, applying beam size correction factors and an additional 
in-flight correction of the absolute response ( 1/1 . 1 in the red and 
1/1.3 in the blue PACS bands) as currently recommended for 
PACS (Poglitsch et al. 2010). The applied beam size correction 
factors are listed in Table [T] 



4. Results and Discussion 

We have detected [OI]63jum and [OIII]52/zm in MIPS J1428, 
and tentatively [O III]52//m in F10214. The spectra are shown in 
Figure [1] (after re-binning to approximately the spectral resolu- 
tion element, which is appropriate for resolved lines). These are 
the first detections of these lines in galaxies at such redshifts, i.e. 
at the epoch of the peak of cosmic star formation. Flux values, 
upper limits and continuum flux densities are given in Table [T] 
Flux uncertainties from the calibration and from uncertainties in 
the continuum definition and line shapes are of the same order, 
and we estimate a total flux uncertainty of 40-50%. 

F10214: we did not detect the targeted lines of [S III]33//m 
and [OIV]26jum, but [O III]52//m (which is observed in the 1st 
order PACS band in pai-aUel with the 2nd order [O IV] data) 
is serendipitously detected at a ~3cr-level. We resolve the ten- 
tative line and measure a FWHM of 300±100km/s (corrected 
for the instrumental profile, which has a FWHM~220km/s at 
170//m), similar to the average CO line width (246±10 km/s, Ao 
et al. 2008). Figure [2] shows the ratio of the limit on [OIV] to 
FIPQ compared to a collection of local template objects from 
the literature (Gracia-Carpio et al., in prep.). The FIR lumi- 
nosity (Lf/R[40-500yum]) has been calculated integrating SED 
data from the literature (Ao et al. 2008 and references therein) 
and our PACS continuum measurements, with a lens magni- 
fication correction factor of 12. The upper limit on [OIV] is 
surprisingly low given that the source is classified as an AGN 
in the optical and mid-IR. Teplitz et al. (2006) report a detec- 
tion of the [NeVI]7.6//m line in their Spitzer-VRSi spectrum. 
Assuming an average [Ne VI]/[0 IV] ratio of 0.4-0.6 in Seyfert 
galaxies (Sturm et al. 2002) and AGN-ULIRGs (Veilleux et al. 
2009) we would have expected a two to three times higher line 
flux in [OIV] than our upper limit. Only a few objects have 
a [Ne VI]/[0 IV] ratio larger than 1 . Given the uncertainties in 
both the Spitzer and Herschel measurements our non-detection 
might still be consistent with a [Ne VI] detection. Figure |2] in 
any case, shows that the upper limit on the [O IV]/FIR ratio in 
F10214 is lower than in Seyfert galaxies, and more in the regime 
of ULIRGs and starbursts. If the differential magnification model 
(see section |2]l is applied, then the upper limit on [O IV] is yet 
lower by a factor of ~3. The low [OIV]/FIR ratio in F10214 
is consistent with L^/r arising predominantly from the starburst 
rather than the AGN, an assumption that is supported by model 



' HIPE is a joint development by the Herschel Science Ground 
Segment Consortium, consisting of ESA, the NASA Herschel Science 
Center, and the HIFI, PACS and SPIRE consortia. 

- FIR(42-122yum)=1.26xl0-"'(2.58x/60 + /100) [Wm^^]. 
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Table 1. Observation details and measured line fluxes 



Source 


line 


OD 


OBSID 


n„,,xn,"„. 


aor duration 
[sec] 


flux" 
[10-"*W/m2] 


continuum flux 
density [mJy] 


corr."^ 


MIPS J 1428 


[O III] 52 


205 


1342187779 


2x7 


5348 


3.7 (0.8) 


117(35) 


1/0.62 


MIPS J1428 


[0 I] 63 


205 


1342187779 


1x7 


2730 


7.8(1.9) 


168 (50) 


1/0.54 


MIPS J1428 


[0 III] 88 


205 


1342187779 


7x7 


18004 


_d 


_d 


- 


F10214 


[0 IV] 26 


179 


1342186812 


10x4 


24827 


<6'= 


330 (100) 


1/0.69 


F10214 


[S III] 33 


185 


1342187021 


6x4 


5864 


<3'= 


378(115) 


1/0.65 


F10214 


[O III] 52 


179 


1342186812 


simult. 


with [O IV] 


0.9 (0.3) 


445 (130) 


1/0.47 



Number of line/range repetitions per nod cycle x number of nod cycles; '' line and continuum fit uncertainty in brackets; an additional 
calibration uncertainty of 30% applies; "^ Correction factor, applied to line fluxes and continuum flux densities, to account for PSF losses of 
the central spatial pixel (see text); "^ The signal in this PACS wavelength range (i.e. >190//m) is contaminated by a 2nd order leak, making 
estimates of upper limits and RMS unreliable; "^ assuming a Gaussian profile with 3-cr peak height and FWHM=300km/s. 
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Fig. 2. The [O IV]/FIR limit in F10214 compared to template ob- 
jects compiled from the literature (Gracia-Carpio et al., in prep.). 
Black symbols: AGN and ULIRGs known to harbour an AGN; 
blue: HII galaxies; green: LINERs. Open symbols and arrows 
are upper line flux limits. 



fits to the FIR SED, as F10214 has strong emission from the FIR 
to the (sub-)mm. 

The [OIIIj/FIR ratio for F10214 is shown in Figure [Sji. It 
is at least an order of magnitude lower than local starbursts 
or AGN. There are not many ULIRG data points available for 
comparison. We have extracted two upper limits on Arp220 
and Mrk231, which are shown in the figure, from ISO-LWS 
data. Initial results from SHINING (Fischer et al. 2010) indi- 
cate a general deficiency of FIR fine-structure lines (including 
[OIII]) in local ULIRGs, similar to the [OIII] deficiency seen 
in F10214. With future data from SHINING we will be able 
to populate this diagram with many more data points for local 
templates covering a large parameter space, in order to put the 
findings reported here into context. 

MIPS J1428: the [OI] and [OIII] Unes (Figure [B are de- 
tected with ~5 cr significance. They are redshifted by 220+100 
km/s and 300+135 km/s, respectively, with respect to z=1.325 
(derived from Ha and CO). Both lines are resolved with a 
FWHM of 750+150km/s each (corrected for the instrumental 
profile). This is considerably broader than CO (380+100 km/s, 
lono et al. 2006), but comparable to Ha (530+160km/s, Borys et 
al. 2006) within the errors. The measured continuum flux densi- 
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Fig. 3. [O I]63yL(m and [O III]52;um fine strength in F10214 
and MIPS J 14284 as function of star formation efficiency. 
Comparison data are from the literature (Gracia-Carpio et al., in 
prep.), exceptfor [OI] in Mrk231 (Fischer et al. 2010). Symbols 
as in Figure |2l 



ties at 120 and 150 fim {Tah.\T^ are low compared to the MIPS 
160 yum flux (430±90mJy, Borys et al. 2006). However, prelim- 
inary processing of PACS and SPIRE photometry of the source 
(S. Oliver and H. Aussel, private communication) is consistent 
with our spectroscopic values. Using the PACS spectrum we de- 
rive Lf/R(40-500//m)~1.3xlO'"'Lo///, where fi is the magnifica- 
tion factor (i.e. L^/r >1.6xlO'^LQforjt/ <8). Contrary to F10214 
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the [OIII]/FIR ratio in MIPS J1428 (Figure|3]l is of the same or- 
der as in local star forming and AGN galaxies. The [OIII]/[OI] 
ratio is very similar to the ratio in typical starburst galaxies (like 
M82). 

For the [O I] spectrum we cannot rule out an underlying ab- 
sorption of ammonia (NH3). In Figure [1] top panel, the dot- 
ted line is a NH3 model spectrum (arbitrarily scaled). The line 
at the bottom of panel (a) is the residual after subtracting the 
line (-Hcontinuum) fit and the NH3 model from the data. Strong 
FIR NH3 features with high column densities have been de- 
tected before in infrared bright galaxies, in the infrared (e.g. in 
Arp220, Gonzalez- Alfonso et al. 2004) and at mm wavelengths 
(e.g. Henkel et al. 2008). On the other hand, the corresponding 
NH3 column density in MIPS J 1428 would be quite high, and the 
object is not known to be heavily obscured at other wavelengths. 

[O I] and [C II] are usually the brightest cooling lines of the 
cool ISM in galaxies. However, ISO data have shown the [CII] 
line in local ULIRGs to be about an order of magnitude lower 
relative to the FIR continuum than in normal and starburst galax- 
ies (e.g. Malhotra et al. 2001, Luhman et al. 2003). The relative 
weakness of [C II] has strong implications on its potential use 
as a star formation tracer for high redshift (z>4) studies (e.g. 
Maiolino et al. 2005 and 2009). Hailey-Dunsheath et al. (2010) 
have observed [C II] with the ZEUS spectrometer at the CSO. 
They found MIPS J 1428 to be a counter example, with a normal 
[C II]/FIR ratio. Our [OI] observation (Fig. [3]) supports this, i.e. 
this source, a high redshift source with the high SFR and L/« of a 
local ULIRG, does not show a deficit in the major PDR cooling 
lines relative to the FIR. 

We plot in Figure |4] the position of MIPS J 1428 in a 
[CII]/[OI] versus ([CII]-h[OI])/FIR diagram. For reference, we 
overplot the PDR model curves from Kaufman et al. (1999), and 
show the location of various samples of template galaxies from 
the literature. This is the first time that a combined [CII] and 
[O I] PDR diagnostic is possible for a source at redshift ~ 1 . The 
PDR model curves are appropriate for clouds in the active re- 
gions of galaxies that are illuminated on all sides (see Contursi 
et al. 2002). Combined HII region and PDR modeling of normal 
and starbursting galaxies have demonstrated that PDRs account 
for more than half of the [CII] emission in these sources (e.g. 
Colbert et al. 1999). Strong shocks could boost the [OI] emis- 
sion relative to [CII]. However, the measured [CII]/[OI] ratio 
of 2.5+1.2 is well above 1, suggesting we can ignore the role of 
shocks in this case. We conclude that local ULIRGs and AGNs 
are not good analogs of MIPS J1428, and that this source appears 
to be more like a normal/starburst galaxy. 

Typical massive galaxies in the distant Universe formed stars 
an order of magnitude more rapidly than in the local Universe. 
Either star formation was significantly more efficient or these 
young galaxies were much more gas rich (Tacconi et al. 2010). 
Figure [3] shows the [O III]/FIR ratio as a function of star for- 
mation efficiency (SEE, expressed in Lfm/Mgas)- Mgas was cal- 
culated from CO measurements (lono et al. 2006, Ao et al. 
2008). For F10214 this assumes a CO conversion factor like 
in local ULIRGs, and that the FIR arises from star formation 
rather than an AGN, an assumption that is supported by the low 
[O IV]/FIR ratio and model fits to the FIR SED (see above). For 
MIPS J 1428 we used a CO conversion factor for normal star 
forming galaxies, motivated by our findings above, which yields 
Mgas = SxlO^^MJfi. The stai- formation rate in MIPS J1428 of 
26OOM0/yr/jU (i.e. SFR > 300Ma/yr for fi < 8), is quite high, 
comparable to local ULIRGs or to SMGs. However, given the 
larger gas reservoir, the star formation efficiency is well in the 
range of normal star forming galaxies. Even with a ULIRG- 
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Fig.4. [CII]/[OI] vs. ([CII]-H[OI])/FIR. PDR models of 
Kaufman et al. (1999) are used. The position of MIPS J 1428 
is indicated, using [C II] from Hailey-Dunsheath et al. (2010). 
Symbols as in Figure |2] sources with L/r >1O'^L0 are marked 
with a cross. Open symbols correspond to [OI] non-detections. 
The Mrk 231 data are from Fischer et al. (2010). 



like conversion factor, LfjR/Mgas would still be lower than in 
most local ULIRGs. Irrespective of the conversion factor we find 
Lfir/Lco = 115Lq/(K km/s pc^) for MIPS J1428. In a recent 
analysis of CO data from z~l-3 normal star forming galaxies 
Genzel et al. (2010) find values of 5 to 100 for this ratio, while 
the average ratio for major mergers is ~160, reaching values up 
to 600. This is consistent with our interpretation of MIPS J 1428 
as a normal star forming galaxy. Our result thus adds further 
support to the recent results from the SINS survey (Tacconi et 
al. 2010) which indicate that with an increased gas reservoir star 
forming galaxies at high redshifts can achieve ULIRG luminosi- 
ties without being major mergers. 

Upcoming SHINING and other Herschel observations will 
enable us to compare with more template objects covering a 
larger parameter space and to search for galaxy evolution eff'ects 
in the ISM of infrared bright galaxies. 
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